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1. Introduction 

The space-time evolution of a quantum mechanical system is governed by its 
Schrodinger equation and in its most complete form it is encoded in the propagator. 
The propagator defines the probability amplitude for a particle to move from one 
point of the space to another in a given time. Similar to propagators in relativistic 
field theories, it provides a global picture of the causal structure of a quantum system 
which goes beyond the information contained in a single wave function. Moreover, it 
plays the essential role in solving the probability related Cauchy problem of Quantum 
Mechanics (QM). 

The vast literature on QM propagators, summarized e.g. in [1], lists mainly 
explicit expressions of propagators for stationary Schrodinger equations in one space 
dimension which are reducible to hypergeometric differential equations (DEs) or their 
confluent forms. This is in strong contrast to the available methods of supersymmetric 
Quantum Mechanics (SUSY QM) [51 [3] which usually lead to much broader classes of 
exactly solvable Schrodinger equations with solutions in particular expressible in terms 
of linear combinations of hypergeometric functions (see e.g. [J). From the structure 
of the SUSY induced relations between superpartner Hamiltonians it is clear that via 
the corresponding Schrodinger equations these relations should extend to relations 
between the associated propagators. The main goal of the present paper is to analyze 
these relations between SUSY partner propagators and to reshape them into user 
friendly general recipes for the construction of new propagator classes. Our interest in 
such new exact propagators is less motivated by their mere existence or their technical 
subtleties, but rather it is in their applicability to concrete physical problems (see e.g. 
[5]) and here especially to models with well-tailored new properties and to setups 
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which up to now were not related to SUSY techniques at aU (see e.g. [B]). In order to 
derive the corresponding technical tools we concentrate in this paper on the general 
approach which allows to establish the link between the propagators of any two SUSY 
partner Hamiltonians. 

Our main idea is the following. For a given Hamiltonian ho, whose Schrodinger 
equation as a second order differential equation is exactly solvable, the general 
solutions of the Schrodinger equation of its SUSY partner Hamiltonian Hn can be 
explicitly constructed. Knowing the relations between the solutions of the Schrodinger 
equations for ho and one can expect to derive transformation operators relating the 
associated propagators. At a first glance, this approach seems to be technically trivial. 
But when one recalls that the calculation of a closed expression for a propagator 
is usual more difficult than the derivation of a corresponding wave function, one 
may expect the problem to be connected with rather nontrivial technical subtleties. 
Moreover, if one wants to establish the explicit link between the propagators one may 
imagine that sometimes the problem may even become unsolvable. 

Probably the first indication that the problem may have a solution was given by 
Jauslin [7] who constructed a general integral transformation scheme simultaneously 
for Schrodinger equations and for heat equations, but who didn't provide a discussion 
of convergency and divergency of the derived expressions. For the sake of convergency 
he applied his technique to the heat-equation-type Fokker-Planck equation only. In 
general, this result may be extended via Wick rotation to propagators for Schrodinger 
equations of a free particle and a particle moving through transparent potentials. 
But the question of convergency and with it the question of solvability remains to be 
clarified. Another indication that the problem may be solvable has been provided by 
Refs [5J [5] where a similar model has been analyzed at the level of Green functions 
of stationary Schrodinger equations. In the present paper we carefully analyze the 
problem along the ideas first announced in [lO] . 

The material is organized as follows. In the next section [21 we briefly recall some 
main facts and notations from SUSY QM necessary for our subsequent analysis and 
we sketch the definition of the propagator. Section [3] is devoted to the basic tool 
of our approach — the interrelation between propagators for models with partner 
Hamiltonians which are linked by first-order SUSY transformations. In Section |4] 
we generalize these results to polynomial supersymmetry. The general technique 
developed in these sections is afterwards used in section [5] to derive the propagators 
for SUSY partner models of "a particle in a box" , of the Harmonic oscillator and of 
the free particle (i.e. for models with transparent potentials). Section [5] concludes the 
paper. 

2. Preliminaries 

The subject of our analysis will be a polynomial generalizatiorQ of the simplest two- 
component system of Witten's non-relativistic supersymmetric quantum mechanics 
fi2\ [21 [3] described by the Schrodinger equation 

{ildt- H)^{x,t) ^0 xe{a,b) (1) 



) For a recent review see e.g. [Tlj . 
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where H is a diagonal super-Hamiltonian consisting of the two super-partners ho and 
/iTV as components 

f ho \ 

H^i ho.N^-dl + Vo,N{x). (2) 

V hN J 

The interval (a, b) may be both finite or infinite. For simplicity we restrict our 
consideration to a stationary setup with ho and h n not explicitly depending on time so 
that the evolution equation ^ reduces via standard substitution "^{x, t) — ^{x)e~^^*' 
and properly chosen boundary conditions to the spectral problem 

H'^ = E^. (3) 

The time evolution of the system may be described in terms of a corresponding 
propagator. 

Furthermore, we assume that the partner Hamiltonians ho and h^ are intertwined 
by an A'^th-order differential operator L with the following properties: 

1. intertwining relations 

Lho = h^L hoL^ = i+Zi^v (4) 

2. factorization rule 

L+L = PN{ho) LL+ = PN{hN) 
Pn{x) = {x - ao) ■ . - {x - ajv-i) 

Im(ai) = ^ ak^i i,k = 0, . . . , N - 1 . (5) 

Here, the adjoint operation is understood in the sense of Laplace (i.e. as formally 
adjoint with the property 9+ = —dx, (AB)^ = B^A~^ and i+ = — i) and the roots 
ai of the polynomial play the role of factorization constants. For simplicity we 
assume that the polynomial P^ has only simple roots. The intertwining relations 
together with the factorization rule can be represented in terms of the polynomial 
super-algebra [T^fTTl lH] 

Q2^(Q+)2 = [Q,H]^[Q+,H]^0 QQ+ + Q+Q ^ Pn{H) (6) 
with nilpotent super-charges 

L J ^ \0 

Although the component Hamiltonians ho and hi^ enter the super-Hamiltonian 
([2|) in an algebraically symmetric way, we consider ho as given Hamiltonian with known 
spectral properties and /iat as derived Hamiltonian with still undefined spectrum. 
More precisely we assume Vo{x) to be real- valued, continuous and bounded from 
belo-wol so that the differential expression ho — —d^ + Vo{x) defines a Sturm-Liouville 
operator which is symmetric with respect to the usual C'^(a,b) inner product. The 
corresponding functions "0 S C'^{a,b) are additionally assumed sufficiently smootlF)!. 
e.g. ip € C^(a, 6), over the interval (a, 6) C R. Moreover, we assume Dirichlet 
boundary conditions (BCs) for the bound state eigenfunctions of ho, i.e a domain 
Viho) ■= {iP : iP e C^{a,b) n C2(a,5), ?A(a) = iP{b) = 0} (see, e.g. [H]) and the 



^) We assume the potential Vo{x) bounded from below and short ranged, because we restrict our 
attention here to models with a finite number of bound states only. For spectral problems on the 
half-line a repulsive singularity at the origin not stronger than £(£ + £ = 0, 1, ... is possible. 

^) As usual, C^{a, b) denotes the space of twice continuously differentiable functions. 
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operator itself being at least essentially self-adjoint (with a closure that we denote 
by the same symbol /iq)- As usual, eigenfunctions which correspond to the continuous 
spectrum of ho are supposed to have an oscillating asymptotic behavior at spatial 
infinity. Concentrating on physically relevant cases we restrict our attention to the 
following three types of setups: 

(i) : The interval (a, 6) is finite |a|, |6| < cx) so that ho has a non-degenerate purely 
discrete spectrum (see e.g. [15]). 

(ii) : For spectral problems on the half-line {a = 0,b — oo) we consider so called 
scattering (or short-ranged) potentials which decrease at infinity faster than any finite 
power of X and have a continuous spectrum filling the positive semi-axis and a finite 
number of discrete levels; the whole spectrum is non-degenerate. 

(iii) : For spectral problems on the whole real line, (a = —oo,b = oo), we consider 
confining as well as scattering potentials. Confining potentials produce purely discrete 
non-degenerate spectra (see e.g. [T^), whereas scattering potentials lead to two- fold 
degenerate continuous spectra filling the whole real line and to a finite number of 
non-degenerate discrete levels (see e.g. |T5]). 

Everywhere in the text we choose real- valued solutions of the differential equation 
{ho — E)tjj = 0. This is always possible since Vo{x) is supposed to be a real-valued 
function and we always restrict ourselves to real values of the parameter E. 

The intertwiner L is completely described by a set of N transformation functions 
Un (x) , which may be both "physical" and "unphysical'lf}! solutions of the stationary 
Schrodinger equation with ho as Hamiltonian: 

hoUn = anUn n = 0, . . . , TV - 1 . 

In our case of a polynomial Pn{x) with simple roots (i.e. ai 7^ ak) the action of the 
intertwiner L on a function / is given by the Crum-Krein formula |17[ I18| 



^ W{uo,ui,...,UN-iJ) j,^^ 



W{uo,ui, . . .,UN-i) 
with W denoting the Wronskian 

W = W{uo, Ml, ... , Mjv_i) = 



Mo Ml ... Utv- 



U 



N-1 



{N-1) {N~l) {N-1 
Mg ... MjY_i 



(8) 



It links the solutions (j> and tp of the Schrodinger equations with h^ and ho as 
Hamiltonians by the relation (jj = Lip. In particular, for iV = 1 (a first-order 
transformation) ([7]) reads 

(/) = Lf = {-uox/uo + d.x)f = (9) 

Mo 

with W{uo) = Uo- Furthermore, the determinant structure ([7]) of the operator L leads 
to the immediate implication that it has a nontrivial kernel space KerL spanned by 
the set of transformation functions m„: 

KeiL — span{Mo, . . . , wjv-i} dim(Ker_L) — dim(KerL^) = N . 



^) By "physical" solutions we mean solutions belonging to the domain 'D(ho) = {?/) : -)/) g 
C^{a,b) n C^{a,b), ip{o-) = = 0}. All solutions i/j ^{ho) corresponding to a spectral 

parameter E outside the continuous spectrum are interpreted as "unphysical" . In the present paper, 
eigenfunctions corresponding to the continuous spectrum are not used as transformation functions 
and, although physically meaningful, we excluded them from our classification scheme of "physical" 
and "unphysical" solutions. 
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The solutions w„ of the equation hpfVn = anVn are elements of the kernel space of the 
adjoint operator and can be obtained ai 

Wn{uo,Ui, . . .,UN~l) 



KerL^ = spanjwo, ■ • • , vn-i} 



n = 0,...,7V-l (10) 



where Wn denotes the Wronskian built as a determinant of the (A^— 1) x (A^— 1) matrix 
with the u„— related column omitted Wn — W{uo, wi, . . . , . . . , uat-i). The 

potential Vjy of the Hamiltonian h]y can be expressed as \T7\ [T5] 

Vn = Vo - 2[\iiW{uo,ui, . . . ,UN^i)]" • (H) 

In general, it is not excluded that the transformation operator L may move a 
solution of ho out of Ho = C^{a,b) thus transforming a physical solution of ho into 
an unphysical solution of h^- Moreover the inverse scenario is also possible, i.e. L 
may transform an unphysical solution of ho into a physical solution of /iat. In such 
cases the point spectrum of ft. at will differ from that of ho- Subsequently, we will 
concentrate on mild transformations L, which leave most of the original spectrum 
invariant with exception of a finite number of spectral points — a characteristic feature 
of differential intertwining operators L leaving the boundary behavior of the solutions 
of the Schrodinger equation unchanged. Recently in [19j a conjecture has been proven, 
which was originally formulated in [SUj and which states that any iVth-order mild 
differential transformation L can be constructed as a superposition from only first- 
and second-order mild transformations. In this case it is possible to show [3] that 
for problems formulated over the whole R (for infinite values of a and b) Vn behaves 
asymptotically like Vq. Therefore the operator hN is also essentially self- adjoint and 
"lives" in the same Hilbert space Ti as ho- Moreover, since the point spectrum of the 
self-adjoint Sturm-Liouville problems that we consider is non-degenerate there is no 
way to create a new discrete level at the position of an already existing discrete level 
and by this means to increase the geometric multiplicity of that levePi 

According to [18] the necessary condition for an A^th-order transformation to 
produce an essentially self-adjoint operator hN i J^I 

{E~ao).-.{E-aN~i)>Q yEespcc{ho). (12) 

This criterion ensures the mildness of the transformation L leading only to changes 
in maximally N spectral points (of the point spectrum). Specifically, the spectrum 
of hN may contain p points more and q points less than spec{ho), where necessarily 
p + q^N. 

In the present paper we will consider the following possibilities. 



^) See also l lA.lll l in [Appendix A| 

®) This is in contrast to SUSY intertwined non-self-adjoint operators ho, h2 for which a second- 
order intertwiner L can map two distinct discrete levels of ho into a second-order Jordan block of h2, 
i.e. an eigenvalue of geometric multiplicity one and algebraic multiplicity two (for the details see Ref 

The basic idea can be understood as signature preservation of Hilbert space metrics, i.e. V-E 
belonging to the point spectrum of ho the eigonfunctions with = {ipEti's) > should 

map into corresponding eigenfunctions (j>E of hjv with = ||I/V-b|P = {'4'E , LrpE) = 

{tp E , Pn {ho)4' e) > what via ^ implies ||12|I . Here, the equality takes place for those tpE for 
which the point E does not belong to the spectrum of /ijvi i-s. if Pn{E) = and <f>E = Liji^ = 0. A 
more detailed analysis of corresponding sufficient conditions for spectral problems on the whole real 
line in case of scattering potentials and of confining potentials is given in |20| . 
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• The spectrum of /ig is a subset of the spectrum of /ijv- A new energy level may 
be created in the spectrum of if and only if the corresponding transformation 
function u{x) is such that l/u{x) G I?(/ijv)i i-e-, in particular, that l/u(x) is 

integrable and satisfies the Dirichlet BCs. 

• The spectrum of is a subset of the spectrum of ho- An energy level may be 
removed from the spectrum of ho if and only if the corresponding transformation 
function u{x) coincides with the /io~eigenf unction of this level, i.e. when u{x) 
satisfies the Dirichlet BCs. 

• The spectrum of Hq coincides with the spectrum oi h^. In this case non of 
the transformation functions ui{x) nor l/ui(x) should be physical, i.e. satisfy 
Dirichlet BCs on both ends of the interval (a, b). This property should be fulfilled 
for all transformation functions ui{x) from which the transformation operator is 
built. 

In all cases we assume that the transformation functions {ui{x)}f£i are linearly 
independent one from the other and their Wronskian W{ui^ . . . , does not vanish 
Va; G (a, b). 

In the remainder of this section, we briefly recall some basic properties of quantum 
mechanical propagators. Once, the super-Hamiltonian ([2]) is diagonal it suffices to 
restrict to one-component (scalar) propagators. 

It is well known that the propagator K{x,y,t) of a non-stationary Schrodinger 
equation contains the complete information about the space-time behavior of the wave 
function "^(xjt) evolving from an initial configuration 5'(a;,0) — '^o{x) 

*(x,t) = / K{x,y,t)-^iy,0)dy 



solving in this way the probability related Cauchy problem of QM. The propagator 
(integration kernel) satisfies a differential equation with Dirac delta function as initial 
condition 

[idt~h{x)]K{x,y;t)^0 K{x,y;0) ^ S{x ~ y) . (13) 
For non-dissipative systems, like in our case, the propagator K{x,y;t) can be 
interpreted as coordinate representation of the unitary evolution operator U{t): 
K{x,y]t) = {x\U{t)\y), where unitarity implies the symmetry K*{x,y; —t) = 
K{y,x;t). 

Subsequently, we will mainly work with a spectral decomposition of propagators 
in terms of complete basis sets of eigenfunctions 

K{x,y;t)^Y.^"^^'*)'('n{y)+ dkMx.tWM (14) 

n=0 

where summation over the point spectrum and integration over the continuous 
(essential) spectrum are understood. By SUSY-transformations we will only induce 
changes in the point spectrum of so that, for simplicity, we will work with 
decompositions over discrete sets of basis functions (corresponding to point spectra) 
keeping in mind that extensions to the continuous spectrum are straight forward. 

It is clear that the defining equation for the propagator of the non-stationary 
Schrodinger equation with super-Hamiltonian ^ can be trivially decomposed as 

Ko(x,y;t) 

[ildt- H]K{x,y;t) =0 K = 



KN{x,y;t) 



Kix,y;0) = ld{x-y). 
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3. Propagators related by first-order intertwiners 

In this section we study the structure of propagators interrelated by first-order SUSY 
transformations. Such SUSY transformations are generated from a single function 
u{x) and they act as basic building blocks in chain representations of higher-order 
intertwiners. Although any iVth-order intertwiner may be represented as a chain of 
first order intertwiners [3] , one has to distinguish between chains which are completely 
reducible within a given Hilbert space Ti. |22j and chains which are partially or 
completely irreducible in Ti.. Complete reducibility means that apart from Hq and 
/ijv also all intermediate first-order SUSY related Hamiltonians /i^, k 0,N are self- 
adjoint or essentially self-adjoint in the same Hilbert space Ti.. In case of irreducible 
chain^ several or all intermediate Hamiltonians are non-self-adjoint in Ti.. Below 
both chain types will play a role. We note that chain representations lead to extremely 
simplified transformation rules for higher-order intertwined propagators and allow for 
very efficient calculation techniques (see subsection 14.21 below) . 
According to ([7]) a first-order intertwiner has the form 

Lx — ~u'{x)/u(x) + dx fiQU — au . 

Due to _E — a > (see and depending on the concrete form of the function 

u{x) the intertwiner L may result in the following three types of relations between the 
spectra of the Hamiltonians /iq and hi (see, e.g., [2]): 

(i) for a — Eq and u — the ground state level Eq of ho is removed from the 
spectrum of hi, 

(ii) hi has a new and deeper ground state level E-i — a < Eq than ho, 

(iii) the spectra of hi and ho completely coincide (a < Eo). 

In order to create a potential Vi (x) which is nonsingular on the whole interval (a, 5) 3 x 
the function u{x) should be nodeless inside this interval. This property is evidently 
fulfilled for type (i) relations since the function u{x) coincides in this case with the 
ground state eigenfunction u(x) — ^o{x). In the cases (ii) and (iii) the nodelessness 
should be ensured by an appropriate choice of u{x), a choice which is always possible 
because of the 'oscillation' theorem (see e.g. [M:). In case (iii) it implies a < Eo- 
Introducing the Green function 

Go{x,y;E) = E -E ^ ' 

ni—O ^ 

of the stationary /ig^Schrodinger equation at fixed energy Epil (see, e.g. [33]) and its 
"regularized" version 



Goiz,y,Eo)=^ ^Jun^ 



Go{z,y,E) 



Eo-E 



m— 1 

the corresponding structural relations for the propagators can be summarized in the 
following 

Theorem 1. The propagators Ki{x,y;t) and Ko{x,y;t) of non- stationary 
Schrodinger equations with SUSY intertwined Hamiltonians hi and ho are interrelated 
with each other and with the Green functions Go{x,y; E) and Go{z,y, Eo) in the 



For a careful analysis of different kinds of irreducible transformations we refer to I19| . 
^) We assume that /iq has a purely discrete spectrum. As it was already stated in the Preliminaries, 
a generalization to the continuous spectrum is straight forward. 
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following way: 
Type (i) relation 

Ki{x,y,t) = L^Ly Ko{x,z,t)Go{z,y,Eo)dz . (16) 

J a 

Type (ii) relation 

-iat 



Ki{x,y,t) = L^Ly KQ{x,z,t)GQ{z,y,a)dz + (l)^i{x)4>-i{y)e . (17) 

J a 

Type (iii) relation 

Ki{x,y,t) = L^Ly Ko{x,z,t)Go{z,y,a)dz . (18) 

J a 

Proof. We start from the type (ii) relation and represent the propagator Ki{x, y; t) in 
terms of the basis functions (t)m{x, t) of the Hamiltonian hi (cf. p4|l ). We note that the 
exphcit time-independence of hi imphes a factorization (j)„i{x,t) — (/)m(a;)exp (—lEmt) 
with (j)m{x) purely real- valued. Expressing </),„ in terms of the corresponding wave 
functions of the Hamiltonian /iq, (pm = NmLtp„i, with N,n = {E — a)~^/^ a 
normalization constant (see, e.g. [2])) we arrive at 



Ki{x,y,t)= (t>,n{x)(t>,n{y)e 

m— — 1 



= L^Ly — e -|-(/)_i(x)0_i(?/)e . 

m=o " 

Modulo a normalization factor N, the wave function 0_i of the new ground state is 
proportional to the inverse power of the transformation function u{x), 0_i = N /u{x). 
It remains to express the time-dependent phase factor in terms of the propagator. 
This can be easily done using the evident property of the bound state solutions of the 
Schrodinger equation 

(•b 

Ko{x, z, t)^^{z)dz = ^/.„(a:)e-'^"* (19) 
so that the previous equation reads 

Ki{x,y,t) ^ L^Ly r Ko{x,z,t) "^"^^^^"^^^ dz -f </)_i(x)</>-i(j/)e-^"*. (20) 



The sum in this relation can be identified as the Green function (fT5| . Due to 
Em — a > 0, Vi?„i G spec(/io) this Green function is regular V£'m and the proof 
for type (ii) transformations is complete. 

The proof for type (i) and (iii) transformations follows the same scheme. The 
formally regularized Green function GQ{z,y, Eq) in (i) results from the fact that 
the ground state with energy Eq is not present in the spectrum of hi so that a 
sum X]m>o appears and the ground state contribution has to be subtracted from 
Go{x, y,a = Eq). In case of a type (iii) transformation a sum X]m=o over the complete 
set of eigenfunctions appears in (PO)) and no new state occurs. □ 
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We conclude this section by reshaping relation for the propagator of a system 
with removed original ground state, i.e. of a type (i) transformed system. The 
corresponding result can be formulated as 

Theorem 2. For transformations with u{x) — tpoix) the propagator Ki{x, y; t) of the 
resulting system can be represented as 

1 [y 1 

Ki{x,y;t)^ TT^a; / Ko{x, z-t)u{z)dz ^ ——L^ Ko{x, z;t)u{z)dz . (21) 

u(y) Ja u{y) Jy 

First of all we recall that ipo {x) being the ground state function of kg satisfies the 
zero boundary conditions. To facilitate the proof of Theorem [5] we need the following 
two lemmas. 

V'o(^)'0o(;/) ' 



Lemma 1. Ly lim ( Go(z, y, E) ~ ^"^ ^^"^^^ ) = Jini LyGo{z, y, E) . 



Proof. This result follows from the explicit representation of Go{z,y, E) in terms of 
basis functions. On the one hand, it holds 

r r fr^ r M2)My)\ ^ 'ilJn{z)Ly^„{y) 

Ly ^lun^_ [Goiz, y, E) - -^^) = g ^„ _ ^„ 

whereas on the other hand the kernel property (annihilation) of the ground state 
Lt/io = gives 



lim {LyGo{z^y,E))^ lim y VW/^V^ y 

\n— 1 / n— 1 



1pn{z)Lylljniy) ( _ ^n{z)Ly1pn{y) 

□ 



Lemma 2. Let fi{x,E) and fr{x,E) satisfy the Schrodinger eguation 

hofix,E):=-f"{x,E)+Vo{x)f{x,E)^Ef{x.E) x e {a,b) (22) 

and boundary conditions 

fi{a,E)=0 fr{b,E) = 0. (23) 

Let also E = Eq be the ground state level of ho with ipo^x) as the ground state function 
(we assume that ho has at least one discrete level) then 

fi{x,E)Lyfr{y,E) ^ _ M^)Iy^oiz)dz 
eTe„ WifrJi) My)ll^oi^)dz 

fr{x,E)Lyfi{y,E) ^ Mx)l!i^oiz)dz 

E^o Wifrjl) ,Po{y) l!^ ^oiz)dz 

where Ly — —u'{y)/u{y) + dy with u{y) = ipaiy). 

Proof. First we note that according to Q Lyfr{y,E) = W{u, fr)/u. Next, since 
both u = ^0 and fr satisfy the same Schrodinger equation it holds W'{u,fr) = 
{Eq — E)ufr and hence 

Win, fr) = {E- Eo) [ u{z)fr{z, E)dz (26) 
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where we have used the property W{u, fr)y=b = which follows from the BCs for u 
and fr- Via we find 

E — Eq '■^ 



and hence 



LyU{y,E)^ I uiz)fr{z,E)dz (27) 



Lyfriy,E) ^ E-Eq fyU{z)U{z,E)dz 

WifrJi) Mb,E) mb,E)u{y) ^ > 

where it has been used that the Wronskian W{fr,fi) — fr{x, E)fi{x, E) — 
fi{x, E)f^{x, E) is x-independent and can be calculated at x — b where fr{b,E) = 0. 
Since the spectrum of hg is non-degenerate, the ground state function is unique up to 
an arbitrary constant factor and, hence, u(x) = tpolx), fr{x,EQ) and /;(a;,_Bo) have 
to be proportional to each other 

frAx,Eo)^CrMx) (29) 



and for E Eq only the first fraction in ([28|) remains undetermined. The I'Hospital 
rule gives for this limit 

lim l--^ = (30) 

E^EoMb,E) fi(h,Eo) 

where the dot denotes the derivative with respect to E. Making use of ([50]) and 

Mb,Eo)f[{b,Eo)= f?{z,Eo)dz (31) 

J a 

(which we prove below) relation (f^5|l yields 

fi{x,E)LyU{y,E) _ f[{b,Eo) uiz)Mz, E,)dz 

A WifrJO ^'^"'''"^^(^i^o) u{y)f^f?iz,Eo)dz ^ ^ 

and via it leads to the result The proof of (|25p follows the same lines with 

evident changes. 

Finally, it remains to derive equation ([?!]) . This is easily accomplished by 
multiplying Schrodinger equation ([22]) for / = fi{x,E) by fi{x,E), its derivative 
with respect to E by / = fi{x,E), and integrating their difference over the interval 
(a, 6). The intermediate result 

rb 

f?{x,E)dx 

= f[{a,E)fi{a,E) ~ fl{a,E)fi{a,E) ~ fl{b,E)Mb,E) + fl{b,E)fi{b,E) (33) 

reduces to (|31l) via BC (1^51) and its derivative with respect to E (what cancels the 
first two terms) and the limit E = Eq, its implication P5|) and the BC for u{x). □ 
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Proof of theorem For the Green function G{x,y, Eq) in (|16p we use the 

standard representation in terms of two linearly independent solutions /; ^ of the 
/iQ— Schrodinger equation introduced in Lemma (see, e.g. [TS]): 

G(x, y, E) = [fi{x, E)fr{v, E)Q{y ~ x) + fi{y, E)fr{x, E)Q{x - y)] /W{fr, fi) (34) 



where 8 denotes, as usual, the Heavyside step function. Then relation (1161) takes the 
form 



Ki{x,y,t) 

i-b 



LxLy 



Kn(x,z,t) lim 

E^Ea 



WifrJl) ~' ' WifrJl) 

where the step functions can be resolved to give 



Eo-E 



dz 



ry 

Ki{x,y,t) ^ LxLy / KQ{x,z,t) lim 

Ja E^Eo 



fi{z)fr{v) _ iJJoiz)i;o{y) 
Wifrji) Eo-E 
fi{y)fr{z) ipo{z)il)Q{y)' 



dz 



dz . 



+ L^L,, / Kr)(x,z,t) lim 

^4 ^ ' E^Eo[W{frJl) Eo~E 

The second argument of the functions fi^r has been omitted for notational simplicity. 
Explicitly acting with the differential operator Ly on the integrals with variable 
2/— boundary yields 



Ki{x,y,t)^Lx Ko{x,z,t)Ly lim 



E-*Eo 



Lx 



Kn{x,z,t)Ly lim 

E^Eo 



fi{z)fr{y) _ ^po{z)■^|;o{y) 
WifrJi) Eo-E 

fi{y)fr{z) ^o{z)i}a{y)' 



dz 



dz 



(35) 



W{frJl) Eo-E 

whereas via Lemma [1] the intertwiner Ly and the limit Yitcie^Eo can be interchanged 
to give 

Ki{x,y,t) = 



(36) 



Lx 



Ko{x,z,t) lim 



fi{z)Lyfr{y) 



E^Eo W{fr,fl) 

Application of Lemma [5] leads to 



dz + Ko{x,z,t) lim 



fr{z)Lyfi{y) 



E^Eo WifrJl) 



dz 



'ipoiq)dq Ko{x,z,t)%l)o{z)dz 



i'o{^)dq / Ko{x,z,t)iljo{z)dz 



(37) 



which we further reshape by expressing the integral with respect to q over the interval 
(y, h) by the difference of two integrals over the intervals (a, h) and (a, y). Substitution 
ipo — u in the first term results in 



Ki{x,y,t) = r^Lx 



Ko{x, z, t)u{z)dz 
rv 



i}}o{q)dq Lx I Ko{x, z,t)^o{z)dz . 



(38) 



My) la i'oii)dq 

The very last integral is nothing but the ground state stationary wave function 
^oix,t) = w(a;)exp (—iEot). Therefore, since Lxu{x) = 0, we obtain the first equality 
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in (j2ip . The second equality results from applying a similar transformation to the 
second term in (|57|) . □ 
The following remarks are in order. First we have to note that the integral 
representation (|2ip is only valid in case of first-order SUSY transformations which 
remove the ground state level. If one wants to create a level in a problem on the 
whole real line one has to use a transformation function u{x) which diverges for 
X ±00 ensuring in this way the normalizability and Dirichlet BCs of the new 
ground state wave function cx l/u{x). An attempt to calculate the propagator 

Ki via (HH) would usually lead to a divergent integral. The correct approach is to 
use (|17p in this case. Jauslin [7] using a different procedure obtained the same result 
((2T|) both for removing and creating a level, but he completely ignored questions of 
convergence or divergence of the corresponding integrals. In concrete calculations he 
avoided divergent integrals by considering the heat equation only. 



4. Higher order transformations 

Addition of new levels 

Let us consider an iVth-order [N = 2, 3, . . .) polynomial supersymmetry corresponding 
to the appearance of N additional levels in the spectrum of hif compared to the 
spectrum of ho. In this case new levels may appear both below the ground state energy 
of ho (reducible supersymmetry) and between any two neighbor levels of ho (irreducible 
supersymmetry, see e.g. |20j). The propagator for the transformed equation can be 
found in the following way. We develop KN{x,y;t) over the complete orthonormal 
set {(pmixjt)} of eigenf unctions of and express all (pm with eigenvalues already 
contained in the spectrum of ho in terms oiipm, i-e. (t>m — NmLipm- The normalization 
constants Nm for transformations fulfilling condition p2)) have the form [3j: 

= [{E - ao){E -ai)...{E^ a^-i)]-^'^ . 

All other eigenfunctions of /ijv which correspond to new levels and which are not 
contained in spec(/io) we keep untouched. This yields 

K^{x,y^t) = L.Ly tll^^^lM ^e-'^-* + X:^n(^)'/'«(2/)e-^""* ■ 

- ao) . . . [Em - aN-i) ^ 

Here we interchanged the derivative operators present in L^ y and the summation. 
This interchange is justified because the propagators are understood not as a usual 
functions but as generalized functions [24 (which, in particular, may be regular, 
i.e. defined with the help of locally integrable functions). It remains to express 
'!/'m(a;)exp {—\Emt) with the help of (fTO]) in terms of Ko{x,z,t), to make use of the 
identity 



Af-l 




and to represent the sum over m in terms of the Green function Go{z,y,an)- As a 
result, one arrives at 

N-i I N-1 ^ \ r°° 

KM{x,y,t) ^ L.j,Ly^^ \ / KQ{x,z,t)GQ{z,y,an)dz 

„=0 V ,=0,j^n "n/J-oo 
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N-l 

+ ^0„(a;)^„(y)e-'""*. (39) 

n=0 

4-2. Removal of levels 

In this section we need information on the intermediate transformation steps of A'^th- 
order SUSY-transformations which goes beyond that presented in section [21 Therefore 
we start with a more detailed description of the transformation operators and solutions 
of the Schrodinger equation at each transformation step. 

Let us consider a chain of N first-order transformations 

tiN < — h-N-i < • • • n-i < n-o 

built from operators ifc+i,fc which intertwine neighbor Hamiltonians hk and hk+i 
as Lk+i,khk = hk+iLk+i,k- We assume all Hamiltonians hk-, k — 0,...,N self- 
adjoint or essentially self-adjoint in the same Hilbert space H so that the SUSY- 
transformation chain itself is completely reducible. Furthermore, we assume that 
at each transformation step the ground state of the corresponding Hamiltonian is 
removed. This means that after TV linear SUSY-transformations the first N states of 
the /iQ— system are removed and the -I- 1st state of ho maps into the ground state 
of /lAT. In our analysis these first -I- 1 states of ho will play a crucial role and we 
denote them by wq,™, n = 0, . . . , N. Furthermore, we use a numbering for the solutions 
Uk^n of the Schrodinger equations of the SUSY-chain Hamiltonians hk, k = 0, . . . , N 
which is "synchronized" with the level numbering of h 0, meaning that a function tik^n 
is related to the spectral parameter En- We have to distinguish between physical 
solutions, which correspond to the existing bound states of hk and which have indices 
n — k, . . . , N , and unphysical auxiliary solutions Uk^n with n = 0, . . . , /c — 1 which we 
construct below. The ground state eigenfunction of a Hamiltonian hk is given by Uk,k 
and for < it is annihilated by the SUSY-intertwiner Lk+i,k 

Lk+i,kUk,k = . 

The bound state functions Uk+i,n, n — k + 1, . . . , N oi hk+i may be obtained by acting 
with the SUSY-intertwiner 

r / , fi T f W{uk,kJ) 

i^k+l.k — ^Uk,k,x/Uk.k + Ox i^k+l.kj — (,4Uj 

Uk,k 

on the corresponding eigenfunctions of hk 

Lk+l,kUk,n ^ Uk+l,n n = k+l,...,N. 

Next, we note that the chain of fc, (fc = 2, . . . , A^) first-order transformations is 
equivalent to a single fcth-order transformation ([7]) generated by the transformation 
functions uo.OjUq.i, . . . ,uo,fc-i- Furthermore, the transformation operators obey the 
composition rules 

Lk+i,kLk,i= Lk+i,i / = 0,...,fc-l fc = l,...,A^-l (41) 

so that, e.g., the second-order transformation operator Lk+2.k intertwines the 
Hamiltonians hk and hk+2 

Lk+2,khk = hk+2Lk+2,k k = 0, . . . ,N - 2 . 

The A'^th-order transformation operator Ln o is then inductively defined as Ln^ — 
Ln,n-iLn-i,o ■ Obviously, it annihilates the A^ lowest states of the original 
Hamiltonian ho, i.e. uofi, ■ ■ ■ ,Mo,iV-i € KerLAr^o- 
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As further ingredient for the derivation of the propagator-mapping we need 
the set of unphysical auxihary functions UAf.n: = 0, ■ ■ ■ , N ~ 1. We construct 
them as upf^n = Lpf QUQ „ , where the functions Mo,n are the unphysical solutions 
of the /iQ— Schrodinger equation at energies En which are linearly independent from 
the eigenf unctions UQ^n- Normalizing wo.n by the condition W^(uo,tu ^o.n) — 1 Sind 
integrating this Wronskian gives 

dy 



and finally 



Ul^niy) 

dy 



In [Appendix A| we show that 



UN.nix) = LN,0U0.n{x) 2 /.A u = 0, . . . , N - 1 . (42) 

J Xo ^ 

Wniuo,...,UN-l) 

UN,n = CAr,„— — ^ (Ad) 

W{uo, . . . ,UN-l) 

CN,n — (En-1 — En){E]y^2 — En) ■ ■ ■ {En+1 — En) n = 0, . . . , iV — 2 

WN-liuo,...,UN-l) W{ua,...,UN-2) 

UN.N^i = — TfTT ^ TJ77 r • (44) 

W{uq,...,un-i) W(uq,...,un-i) 

The ground state function oi is obtained by acting with Ljv.o on the A^th excited 
state of Hq: 

W{uq,q,...,Uo^n) 

un.n = En,oUo,n = 7777 7 ■ (45) 

W(Uo,o,...,Uo.N-l) 

Below we will often use the abbreviations W{x), Wn{x) for the Wronskians (when 
there is no ambiguity in the definition of the transformation functions wo,m in their 
arguments) indicating explicitly only the dependence on the spatial coordinate. 

Within the above framework, the A^ first discrete levels Eq, Ei, . . . , -Eat-i have 
been removed from the spectrum of ho by choosing the ground state functions Uk,k of 
the Hamiltonians as intermediate transformation functions. For such a construction 
the transformed propagator may be calculated according to 

Theorem 3. Let the N first eigenf unctions uo^n = Un = i^n, n — . . . , N — 1 of h^ he 
the SUSY transformation functions. Then the propagators Kn{x, y; t) and Kq{x, y; t) 
of the Schrodinger equations with Hamiltonians h pf and h^ are interrelated as 

N-l 



Kn{x, y- 1) = (-l)^LAr,o,. E WT ^' ^)"»(^)'^^ (^6) 

„=o ^^y> 

= (-l)^-iLiv,o,x (-1)" WT / ^' *)"»(^)^^ • (47) 



Proof. The proof of these relations can be given by induction. We start with (|47l) . For 
Ki{x,y;t) the statement is proven in (I2ip . Assuming that (|47|) holds for KN{x,y;t) 
we verify its validity for fCjv+i(a;, y; t). The corresponding Hamiltonians ft-Tv+i and hj^ 
are intertwined by the linear transformation L^j^i^n so that (I2ip is applicable and 
Kn+i{x, y; t) can be represented as 

1 

KN+i{x,y]t) ^ i-^Ln+i^n.x I KN{x,z;t)uN,Ndz . 

UN,N[y) Jy 
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Replacing K^^ by l|i71) and making use of relations ([33]), pi)) and the composition rule 
(Plj) gives 



(-l)^^-^ifw+i(x-,y;t) 



1 



UN,Niy) 



L 



N + l,0,x 



N-1 



X ^(-1)"C^W dz / d(7MAr^„(z)uAr,Ar(z)ii:o(x,g;t)wo,„(g) . (48) 



The integration region of the double integral is the upper triangle of the rectangle 
y < z,q < b in the {z, g)— plane. We replace this double integral by the difference of 
two double integrals over the whole rectangle and the lower triangle, respectively, 



N-l 



{-lf-'KN+,ix,y;t) = ^(-l)"C^i 

UNMy) ^0 



dzKQ{x,z;t)uQ^n{z) I dquN,n{q)uN,N{q) 



dzKo{x,z;t)uQ^n{z) / dquN,n{q)uN,N{q) 



(49) 



Here, we reshape the two integrals of the type L dquN,n{q)uN.Niq) as follows. First 
we note that UN.n and un,n are solutions of the same Schrodinger equation with 
Hamiltonian hj^ and therefore 

, / N / N W[uN,N{0,UN,n{0] w 
dqUN,n{q)UN,N[q) = Wb,n 



En — E]^ 

un.nLn+i.nuna 
En - En 



(50) 



where Wb^n ■— W [ujv_jv(6), MAr,„(6)] /{En — E^) and where the second equality was 
obtained via (1^0)1 . Applying the general relation (|43p to LN+i.NUN,n = ttTV+i,™ leads 
finally to 



, , . , . ^ Wn{uo, ...,un) 

dqUN,n(q)UN,N[q) = - CN.nUN.N-rr-, T" - Wb,., 

^ ' W{uq,. . . ,Un) 

= : -CN.nUN.N{0^^-Wb,n. 

With ([5T|) as substitution rule the propagator takes the form 

Wn{y) 



(51) 



{-ir-'KN-,,{x,y-t) = / 



1 



N-l 



Ko{x, z]t)uo,n{z)dz 



LN+ifi,x I dzKo{x,z;t) 



un,n{z) 

^ ' n=0 



Y,{-\Tu^^n{z)Wn{z). 



UN,Niy) 

(The terms containing Wb^n exactly cancelled.) The sum 



Af-l 



Sn ■■= ^ (-l)"-uo,„W„(uo, . . . , ujv) 



(52) 



(53) 



n=0 
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in the second term can be calculated explicitly. Comparison with the evident 
determinant identity 



■"0,0 


• ■ uq,n-i 


U0,N 


■"0,0 


■ ■ ■"0,Af-l 


U0,N 


(N~l) 
0,0 


JN-1) 
■ ■ "o,Af-l 


"o,Ar 



uo,o 

"0,0 


■ • ""0,^-1 

■ • ""0,^-1 




Uo,N 


+ 




"0,0 




• ■ ""0,^-1 


Uo,N 


(N-l) 
0,0 


■ ■ "o,Af-l 


"0,JV 




"o,o 


(N-l) 
■ ■ ""0,^-1 


"0,JV 



(54) 



shows that (|53p coincides with the decomposition of the first determinant in the second 
line of over the elements of its first row. Hence, it holds 



S 



N 



(-1) Mo,JvM^(^"o,o, ■ • • , "0,^-1) 



Representing the ground state eigenfunctions un^n in (|52[) via (|45p in terms of 
Wronskian fractions we find that 



( \ V—. 



and, hence, that the second term in ()52p is nothing but the absent n = N summand 
of the sum in the first term. As a result, we arrive at 



N 



KN+iix,y;t) = (-l)^LAr+i,o,.^(-l)' 



Wnjup, . . . ,Mjv) 
W{uo, . ■ . ,UAr) 



Ko{x,z;t)uo^n{z)dz 



what completes the proof of (|47p . The representation follows from (|¥71) and the 
relations 

rb 

Ko{x, z, t)un{z)dz = M„(a;)exp {~\Ent) LNfl,xUn{x) = . 

□ 

We note that in formulas (j46p and (|47|) only one-dimensional integrals are present. 
In this way, they may turn out more convenient for concrete calculations than similar 
equations derived in [7] . 

Furthermore, we note the following. Theorem [3] is proven for the case when the 
lowest discrete levels are removed from the spectrum of ho starting from the ground 
state level. This scenario corresponds to reducible supersymmetry. In order to see 
which of the conditions on the transformation functions Mo,n used for the construction 
of the propagator representations and (|T7)) are indeed necessary conditions one 
may simply insert Kn{x, y; t) directly into Schrodinger equation (fT3|) . It turns out that 
neither the condition of level deletion starting from the ground state nor a deletion of a 
level block without surviving levels inside is used. This means that equations (|46p and 
(|47p hold for any choice of transformation functions provided their Wronskian does 
not vanish inside the interval (a, b), i.e. it holds for reducible as well as for irreducible 
SUSY transformation chains. A necessary but in general not sufficient condition for 
the nodelessness of the Wronskian is inequality p2|) (for further details see (20]). 
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4-3. Strictly isospectral transformations 

Strictly isospectral transformations can be generated with the help of unphysical 
solutions of the Schrodinger equation as transformation functions. In this section, 
we extend theorem [3] to models defined over the whole real line (a, 6) ~ (— cx3,oo) by 
assuming transformation functions which vanish at one of the infinities a; — > — oo or 
X oo and violate the Dirichlet BCs at the opposite infinities (a; ^ cxi or a; ^ —oo). 

In accordance with (|46|) and (|47|) we formulate the corresponding relaxed version 
of theorem [3] as 

Theorem 4. Let the transformation functions Un{x) vanish at only one of the infini- 
ties X — oo or a; ^ oo of the real axis K. Then the propagators KM{x,y]t) and 
Ko(x,y]t) of the Schrodinger equations with h^ and ho as Hamiltonians are related 
as follows: 

for Un{x — oo) — > : 

Km{x, y; t) = (-l)^L, ^ WT / ^)""(^)^^ (^5) 

„=o 

for Un{x oo) —^ : 

for Uk{x — oo) -^0 fc = 0, . . . , Af and Um{x oo) ^0 m = M + 1, . . . , N — 1 



Km{x, y- 1) = i-lf~'L, (^1)"-^^ *>niz)dz (56) 



i^Ar(a;,y;t) = (-l)^L,^(-l)"-^M f Koix, z;t)ukiz)dz 

+ (-l)^-iL. ^ (_1)™]^lM ri^o(a;,z;t)u„(z)dz. (57) 

Proof. We have to verify that the initial condition KQ{x,y,0) = S{x — y) and the 
Schrodinger equations (idt — hQx)KQ(x, y,t) = Q and (\dt — hoy)Ko{x, y,t) — fulfilled 
by the original propagator KQ{x,y,t) map into corresponding relations for the final 
propagator K]\f{x,y,t), i.e. that Kiq{x,y,Q) = 5{x — y), (idt — hNx)KN{x,y,t) = 
and (idt — hNy)Ki^{x,y,t) = are satisfied. We demonstrate the explicit proof for 
the setup with Uk(x — s- — oo) — s- omitting the technically identical considerations for 
the other cases. 

We start by noticing that the intcrtwiner L,j; maps solutions of the Schrodinger 
equation for ho into solutions of the Schrodinger equation for hj^ and, hence, 
(idt ~ hNx)KN{x, y,t) = is automatically satisfied. 

Next, we consider the initial condition Kpf{x,y;0) — S(x — y) which should be 
fulfilled by the r.h.s. of With Ko{x, y, 0) = S{x - y) and /^^ 5{x - z)u^(z)dz = 
9{y — x)un(x) we have from (|57p 

K^{x,y- 0) = (-1)^ ^^(_i)"^^L^o_^[0(y _ a;)«„(a;)] . (58) 
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W{uo, . . .,UN-l) 

we represent the derivatives d™[9{y — a;)u„(x)], to = 0, . . . , as 

m — 1 
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(59) 



(60) 



fe=0 



^i™ '^\y — x) :— 9™ '^6{'y ~ x). Taking into account that 



uq(x) 
u'q{x) 

,{N) 



UN-i{x) 9{y - x)un{x) 
u'N-iix) 0{y ~ x)u'^{x)) 



ull''{x) ... u']^\{x) 0{y-x)ui!^^ 
and making use of the hncarity properties of determinants we reshape (|59p as 

ijv,o,x[^(a; - y)unix)] 

Uo{x) ... UN-l{x) 

Uq{x) ... u'j^_j^{x) -5{x-y)un{x) 



= 



1 

W{x) 



,(N) 



i(^) J2k=() Q 



^-1 r'Nd^-k) 



{y - x)un\x) 



(61) 



Expanding this determinant with regard to the elements of the last column we find 
LN,oAOix - y)u,,{x)] = ^— ^ E Grei"^~''\y - x)ui^Hx) (62) 

' m=l k=0 

where WVm(x) are corresponding minors. For the verification of the relation 
Kn{x, y; 0) = S(x — y) we use its representation 



(63) 



KNix, y; 0)f{x)dx = f{y) 



where f{x) is a sufficiently smooth test function with compact support. With (|58 
and §^ the l.h.s. of §^ reads 



2^2.(-i) — — y 

n=0 m=l fc=0 



dx e-^^ (y ~ X) WNm{xm-^U^\x) ^g^^ 



As next step, we use 9^ ^\y — x) — — ^ ^\x — y) and multiple integration by 
part^lj to remove the derivatives from the 0— functions: 



N-l N m-1 



l.h.s. Of (El - E E E(-i)""'^^" 



The relation 



n— m— 1 fc— 



AT-l 



w^(y) ^ 



W-Ar,„(7/)/(y)u.L''^(y) 
W^(y) 



(65) 



i_^(_l)»W'„«(f)^(_l)A^<5^.^_^, j=0,...,iV-l (66) 



^) For the theory of distributions (generaUzed functions) see, e.g., |241 (in particular vol.1 p. 26). 
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reduces this multiple sum to 



N 



l.h.s. Of ® = i-lf'^^^^^^^ E (-1)'^^^' 



(67) 



and because of VFjvw(y) = W{y) and Efc=o ("1) ^'fc ^ (^f- 4.2.1.3 in 28]) 

the condition is satisfied. 

It remains to prove that the Schrodinger equation (idt — hNy)KN{x,y,t) = is 
fulfilled too. By explicit substitution of equation (|55p we have 

{idt - hNy)KN(x,y,t 

N-l 



l)^i, E(-1)"WT / ^9tKo{x,z;t)ur,{z)dz 



n=0 
N-1 



Wiy) 



VN{y) 



W„{y) 
Wiy) 



K(){x, z; t)un{z)dz 



n=0 



T„/ N un{y) + — — -(u„(?/)dy + u„(y)) 



W{y) 



Ko{x,y;t). (68) 



First, we note that due to relation (j66p and its derivative the last sum vanishes. Taking 
further into account that Wn{y)/W{y) is a solution of the Schrodinger equation for 
/lAT at energy £"„ (cf. pO)) ) and replacing i9tifo(x, z; <) ^ hozKoix, z]t) one reduces 
equation ([68|) to 



(69) 



(i9t - hNy)KM{x,y,t) 



(- 



l)^i, E(-1)"WT r [{h^z~ E,,)K^{x,z-t)]u^{z)dz 



Integrating by parts and making use of [hoz — En)un{z) = 0, the asymptotical behavior 
u{z — oo) — !■ 0, u'{z —oo) — > and relation (j66p one finds that the r.h.s. 
in vanishes and, hence, the Schrodinger equation (ii9j — hNy)Kpf{x,y,t) = is 
fulfilled. □ 



4- 4- General polynomial supersymmetry 

The three different types of transformations considered above may be combined in 
various ways to produce a supersymmetry of more general type. In general, from the 
spectrum of the original Hamiltonian q levels may be removed and p additional 
levels may be added, p + q < N producing in this way the spectral set oi h^- The 
inequality would correspond to SUSY transformation chains between /ig and which 
contain isospectral transformations. For further convenience we split the spectra of 
/lo and ft. AT according to their transformation related contents as 

spec(/io) = {si, (3j,Ek} + speCc(/io) , i^l,...,q; j ^ 1, . . . , N - {p + q + r) 
spec{hN) = {Xi,f3j,Ek} + spec^{hN) , l = l,...,p; j = 1, . . . , N - {p + q + r) (70) 

where the discrete levels and the continuous spectrum spec^{ho) = speCf,(/iAr) 
are not affected by the SYSY transformations. The set of transformation constants 
{oin}n=o ~ {^i^ Pj^lk} correspouds to p new discrete levels A„ £ spec(/iAr) not 
present in spec(ft.o)j 9 levels G spec(/io) not present in spec(/iAr), N — {p + q + r) 
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levels Pj present in both spectra and r constants 7^, not coinciding with any energy 
level of both Hamiltonians, 7^ ^ spec(/io) U spec(/iAr). Transformations induced at 
constants q;„ = (3mln are strictly isospectral. 

Summarizing the previous results the following expression for the propagator 
Kn{x, y, t) can be given 

N-l I N-l ^ \ .00 

KNix,y,t) ^ LxLy^ Yl / Ko{x,z,t)Go{z,y,an)dz 

+ ^ 0;,„ {x)<t>x^ (y)e-'^"* + ^ 0ft. (x)0A. (y)e-'^"* (71) 
where for a„ = e„,/3„ 

V'n(a;)'0n(y) 



00(^,2;, «„) = lim 



Go{z.y,E) 



an - E 

and Gq{z, y, q;„) = Go(z, y, a„) otherwise. 

5. Applications 
5.1. Particle in a box 

As first application we consider a free particle in a box, i.e. the Schrodinger equation 
with Vo{x) = and Dirichlet BCs at the ends of the finite interval (0,1). The 
eigenfunctions of this problem are the well known ■0n-i(a;) ~ \/2sm{mTx), n = 1,2, . . . 
with energies — n^n'^. The corresponding propagator reads [25] 

Kboxoix,y,t) = i ['&^{x,y;t) - i^j{x,y;t)] 



with 

' x-y 



d^{x,y;t) -.= ^3 



TTt\ „, , . „ fx + y Trt\ 



and 'i?3(g|T) denoting the third thcta function [26] . 

As SUSY partner problem we choose a model which we derive by removing the 
ground state level Eq by a li near (one-step) SUSY-mapping with u = = sinTrcc 
as transformation functior l^^)! . This leads to the Schrodinger equation with potential 
^1(3:) = 27r^/ sin^(7rx), i.e. a singular Sturm-Liouvillc problem. The propagator of 
this problem can be represented via ([2T|) as 



1 n 

Kboxi{x,y,t) ^ -— \-d^{x,z;t)-dj{x,z;t)]sm{Trz)dz (72) 

2smTTy Jo 

or after explicit substitution of Lx = — tt cot(7ra;) + dx as 

Kboxi{x,y,t) = ^'^"^(^^^ / [d-(x,z-,t) - d^{x,z-,t)\su\{-nz)dz 
2sm(7ry) Jo 

y 1 



, [i?3 (x, z; i) + (a^, ^; t)] cos(7rz)dz + - [zJg (x, y; t) + i?+ (x, y ; t)] . (73) 
z sm^TTj/j Jo z 

Here after using the property dx0f{x,z,t) = ±dzOf{x,z,t) we integrated in ([72]) by 
parts. 



There exist other types of transformations leading to regular transformed Sturm-Liouville 
problems. But the solutions of the resulting Schrodinger equations will violate the Dirichlet BCs 
[9] so that a special analysis is needed. 
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Here we consider the Hamiltonian 



X 

T 



Using a two-fold transformation with transformation functions 



U2{x) = (a;^ - l)e" 



usix) — x{x'^ — 3)e 



corresponding to the second and third excited state eigenfunctions of ho we obtain a 
perturbed Harmonic oscillator potential [3] 



8x^ 



96x^ 



X 

T 



(74) 



^4 + 3 {x^ + 3)2 

which for large behaves like the original harmonic oscillator potential, but for small 
|a;| shows two shallow minima. For completeness, we note that the transformation 
functions ^2(3;) and 7/3(2;) have nodes, whereas their Wronskian W{u2,U3) = {x* + 
3)e~^ /2 is nodeless so that the corresponding second-order SUSY-transformation itself 
is well defined, but irreducible. 

The propagator for the Schrodinger equation with Hamiltonian h^^'^'> — + 
y (2,3) ^^-j ^gj^ |-,g constructed from the well known propagator 

1 i[(x^ + y^) co3t-2xy] 

KoSc{X,y,t) _ = e "^int 

v47ri smr 

for the /iQ— Schrodinger equation (see e.g. [57]) via relation ([TO]) . The occurring 
integrals Kosdx, z, t)un{z)dz can be represented as derivatives with respect to the 
auxiliary current J of the generating function 



S{J,x,y,t) 



1 



V 47ri sin t 



cxp 



i[{x^ + z^) cost 
ismt 



2xz] z^ 

--r + J^ 

4 



dz 



1 



-exp 



-2it 



1 + erf 



(ij2 smt -f Jx)exp {—it) 



-■\/iexp (— ^) (2 Jsini + i(?/cxp (it) — x)) 



As a result, we obtain 

K^'-'\x,y,t) 
'y{y^-S) \d^S{J) 



aj2 



S{J) 



2y/smt 



y 



d^S{J) dS{J) 



aj3 



dJ 



J=0 



in terms of obvious abbreviations. 

The technique can be straight forwardly generalized to second-order SUSY- 
transformations built on any eigenfunction pair 

Uk{x) ^ pkix)e~''^/^ Uk+i{x) =pA;+i(x)c"^'/'* 

of the Hamiltonian fiQ. The corresponding generalized potentials V^'''''^^\x) (see [3]) 
read 



Qk{x) 



Qk{x) 



X 



(75) 
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where 

Qk{x) -.^pl^iix) ~ Pk{x)pk+2{x) W{uk,Uk+i) = -Qfc(x)exp (-.T^/2) 

are built from the re-scaled Hermite polynomials Pk{x) = 2^^^"^ Hk{x/ \/2). 
With the help of ([^5]) the propagator is obtained as 

i^C^^'^+i) [X, t) = ey'/^L^ [pk {dj) S{J)] - ^ [pk^^ {dj) S{J)]] (76) 

V Qk{y) Qk{y) J j^o 

with pk (dj) denoting the fcth-order differential operator obtained from Pk{x) by 
replacing a;" 

Finally we note that the potentials ([75)1 behave asymptotically like /A for 
— > oo and have k shallow minima at their bottom. 



5.3. Transparent potentials 



Here we apply our propagator calculation method to transparent potentials which are 
SUSY partners of the zero potential in models defined over the whole real axis. We note 
that the propagator for a one-level transparent potential was previously calculated by 
Jauslin for the Fokker-Planck equation. The propagator for a two-level potential 
can be found in [9] where the general form of the propagator for an A'^-level potential 
has been given as a conjecture. Here we will prove this conjecture. 

We construct the iV-level transparent potential from the zero potential with the 
help of the following N solutions of the Schrodinger equation with Voix) =0 (see e.g. 

W) 

U2j-i{x) = cosh{a2j-ix + &2i-i) (77) 
'U2j{x) = sinh(a2j a; -I- 62^ ) 

as transformation functions. Here, a^, 6^ G ] 
denotes the integer part of (A^ + l)/2 . 

The factorization constants aj — —a: 
(point spectrum) Ej = aj < oi Hn = —d^ + Vn{x), whereas the continuous part of 
the spectrum of /ijv fills the whole real axis. The eigenfunctions of the discrete levels 
normalized to unity have the form [2] 

1/2 

W{ui,U2, . . . , U„-l, U„+l, . . . , Un) 



i-l) 

j = l,2,...,[(7V-t-l)/2] 
a,; > fli+i > is assumed and [(A^+l)/2] 



define the positions of the discrete levels 



iPn{x) = 



N 



V n 



W{U1,U2, .. ., un) 



(78) 



The propagator Kq and the Green function Gq for the free particle are well-known [T] 

1 H--y)^ 



Ko{x,y;t) 



(79) 



Go{x,y,E) 



jK\x-y\ 



so that according to 
as 



KN{x,y,t) 



2k 



ImK > 



E ■ 



the propagator of the transformed system can be calculated 




N 

n 



a„ 



exp 



i(a; — zY 
4t 



an\z-y\ — 



dz 
2a„ 



ipn{x)(pn{y)e' 



Kcn{x, y, t) + KdN{x, y, t) . 



(80) 



Exact propagators for SUSY partners 



23 



In the last line we separated contributions from the continuous spectrum, Kcn{x, y, t), 
from contributions from the discrete spectrum, Kdiq{x,y,t). 

The integral in KcNix,y,t) can easily be calculated since the primitive of the 
integrand is well known (see e.g. integral 1.3.3.17 of Ref [IH])- Using the well studied 
convergency conditions of the error-function erfc from [29) we find 



I{a,x,y,t) := 



1 



4a 
a > 

and with it 



47rit 

e'»(^-a)erfc f aVTt 
t > 



exp 



i(x — z)^ 



At 



dz 

a|.-2/|)- 



x--y\ 
2Vit J 



(81) 



N 



N 



71—1 \m— l,m:^n 



I{an,x,y,t) 



(82) 



Introducing the notation 



erf ± (a) := erf av it =F 



I) 



2^/i^ / 

and abbreviating the Wronskian of the x— dependent transformation functions ui, 
UN as W{x) we formulate the final expression for the propagator as 

Theorem 5. The propagator for an N -level transparent potential has the form 



KN(x,y;t) 



1 



VAttU 



e ft 



N 

E 

n=l 



N 

n I 



Wn{x)W^{y) :2 



W{x)Wiy) 



" [erf +(a„) + erf _(a„)] 



(83) 



For the proof of this theorem we need the additional 

Le mm a 3. Let {ai}fLi be a set of N non- coinciding complex numbers a,; ^ 0'j=ii G C 
an Ji^^i n G Z+ . Then the following identity holds: 



N 



N 



E< n 



1 



ai — a," 



(84) 



Proof. Consider the auxiliary function 



(z - Qi)(z - a2) . . . (z - un) 

which is analytic in any finite part of the complex z-plane except for N simple poles 
ai, . . . , ajv- From the residue theorem follows 



N 



N 



^res/(a,) = 



N 

n 



1 



-res f{oo) 



what with the residue at infinity yields 



□ 



^) We use the standard notation Z+ = {0, 1, 2, . . .} for the natural numbers with zero included. 
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Proof. Without loss of generaUty we may set bj = 0. We start with the propagator 
component Kcn{x, y, t) in (|82p related to the continuous spectrum. First we note that 
the function / in (|82p depends only via the difference x — y on the spatial coordinates so 
that the action of dy in Ly can be replaced by dy (— Hence, the composition 
of the two A^th-order transformation operators L^Ly acts as an effective 27Vth-order 
differential operator in x 



LxLy — i?o + Ridx + . . . -f- rt2Aft 
with coefficient functions Rn{x,y). Accordingly, ([5^ takes the form 



(85) 



N 



N 



KcN{x,y;t)^Y.\ n 



1 



[Ro + Rldx 



R2Ndx ]Jn 



(86) 



n—1 \m— l,m^n 

where Jn — I {a^n t). From the cxphcit structure of the first derivatives of Jn 

dJn _ 1 

dx 4 

d^Jn 



'<^.(:^-y)erfc + (a„) + e°"(^-^)erfc_(a„) 



-1 i(^-y) 



dx"^ V47rii 
a3j„ -1 



= e « + </(a„) 



dx^ 2tV^Tdi 



[X — y)e 



-.2 r 



-e-'^"(^-^)erfc + (a„) +e""(^-^)erfc_(a„) 



(87) 



wc find (by induction) the general expression of an even-order derivative 



g2,nj 



afAk„r{x,y,t)e- 



k=0 



+ 



and of an odd-order derivative 



-""(^-«)erfc + (a„) e''"(=^-^)erfc_(a„) 



Q2m-l7- ™-2 

T = 2^ an Akm{x,y,t)e « 

fc=0 



-'^"(^-^)erfc+(a„)-e°"(^-^)erfc_(a„)l . (89) 



It holds 



Akmix,y,t) = 



m = 2l n fc>Z-l 
m = 2«-M n k> l-l. 



Subsequently we use the abbreviation 

Im{an) - T^e"^"* [(-l)"e-'^"(^-«)erfc + (a„) + e'^"(^-«)erfc-(a„) 
4a„ L 



(90) 
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and we will need the explicit form of the a^*^^ and terms in the highest-order 
derivative 



= ... — a: 



The complete propagator component KcN{x,y]t) in (j86p can now be rewritten as 



(91) 



N 



N 



K^N{x,y;t) n 



1 



n— 1 \m—l^m^n 



-Ro^(an) 



AT 



AT 



n 

n— 1 \ m— l.mT^n 



a„_Ri/i(a„) 



TV 



TV 



N 



E n 

n— 1 \ m— l,my^:n 



E n 

n— 1 \ m— l.mT^n 



R2N 



i?2 



1 i(x-B) 



/Airit 



= e « +a„/2(a„) 



'N-1 



Comparison with Lemma [3] shows that due to = only a single term 

containing exp {\{x — yY /At) does not vanish. It is the k — N ~ 1 term in the very 
last sum which with R2N = (—1)^ and Aiq-i^2N{x,y,t) — — (47rit)~-^/^ yields exactly 
the propagator ((79|) of the free particle. All other terms containing exp (i(a; — ?/)^/4t), 
after interchanging the sums X]n=i ■ • ■ ^'^'^ '^k=o • ■ ■' cancel out due to Lemma [3] 
Thus, the above formula reduces to 

N I N 



KcN{x,y;t) 



1 ,(--yf 
zC <1* 



viTrii 



E n 



1 



(92) 



n— 1 \ m— l,m^n 

X [i?o/o(a„) + i?ia„/i(a„) + i?2a^^2(a„) + Rsalhian) + . . . + i?2wa?i^/2Ar(a„)] 
Expressing /m(an) via (|90p in terms of erfc+ and erfc_ one finds contributions 



+ ari?2jv)e°"(^-^) 



erfc + (an) ig^^ ^ 3 

e " (i?o — a„i?i + a i?2 — a-n^s ' 

4a„ 

erfc (a„) ^i,^,^^^ ^ ^^^^ ^ ^2^^ ^ ^3^^ _ 
4a„ 

which can be represented as 

erfc + (a„) ^i,2,^^^ ^ ^^^^ ^ ^^^2 ^ ^^^3 ^ ^ ^ ^^^^2Ar)^-a„(.-,) 
erfc (a„) ^i,.,^^^ ^ ^^^^ ^ ^^^2 + ^3^3 ^ . . . + i?2jvar)e°"(^-^) • 

4(2 

Comparison with l|85p shows that the sums yield simply LajL^exp [±o„(a; — y)]. 
Recalling furthermore that the transformation operators and Ly are given by ([7]) 
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a nd t hat the transformation functions have the form (|77p we arrive after some algebra 



aO 



L.e±-- = (Tir(-l)^+"-^a„ {[ (a^^ - al)^^ (93) 



and, hence, at 

N 



T T ^±a„{x-y) ( ,\n 2 TT ^ 2 2^2 Wn{x)Wn{y) , . 



Substituting ((Ml) into we obtain 



N ( N 

iEe""W n ■ 

n—l \rn—l,rn^n 



V47rii 



A further simphfication can be achieved by recalhng that Q!„ = — a^^ and that the 
parameters a.^ are ordered as ai < a2 < . . . < a^r- Setting (a^ — a^) = — (a^ — a^) 
in the denominator for m — 1, . . . , n — 1 gives an additional sign factor (—1)""^. The 
propagator sum Kiq — Kcn + KdN resulting from continuous and discrete spectral 
components reads then 

Ti— 1 \ m— l.m^^n 

^ ^(^)W^'(^f [erfc + (a„)+crfc_(a„)] 
and via the relation erfc(z) — 1~ erf(z) it leads to the expression in the theorem. □ 



6. Conclusion 



A careful study of propagators for SUSY partner Hamiltonians is given. Assuming 
the partner Hamiltonians linked by polynomial supersymmetry of a general type 
we derived user friendly expressions interrelating the corresponding associated 
propagators. Since the propagators may be also defined in terms of continual integrals 
the results should be useful in exploring new classes of continual integrals. We applied 
our general technique to derive propagators for transparent potentials and for a family 
of SUSY partner potentials of the Harmonic Oscillator. Further applications can be 
expected in time-evolution scenarios for particles placed in corresponding potential 

We note that equation l|93|l is compatible with H78|l . The function exp {—anx) is a solution of the 
initial Schrodinger equation related to one of the factorization constants and it is linearly independent 
from the corresponding factorization solution. Therefore L^jexp (— a„i) up to a constant factor is one 
of the bound state functions of /i]v given in jTSj. 



Exact propagators for SUSY partners 



27 



fields, in stochastic processes described by Focker-Planck equations (see e.g. [7^), in 
statistical physics and possibly in dynamo processes of magnetohydrodynamics |30| . 
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Appendix A. 

Here we derive a representation of the unphysical auxiliary solutions uat,™, n = 
0, . . . , N — 1 of the /ijv— Schrodinger equation in terms of Wronskian fractions. We 
start by recalling that the functions ujv.n are defined by acting with an A^th-order 
(polynomial) SUSY-transformation operator Lj^^ on the unphysical solutions wo.n of 
the Schrodinger equation with Hamiltonian ho at energies E„ = an- The functions 
are linearly independent from their physical counterparts uq „ and normalized as 
) = 1. The operator Ljv,o itself is constructed from the wo,„, n = 
0, . . . , N — 1 as transformation functions. 

Below we show by induction that the unphysical solutions defined in (|42|) as 

dy 

have the following representation in terms of Wronskian fractions 

^ VF"„(uo,o, • ■ ■ ,"o,Ar-i) / . 

UN,n = L-N.nTTTT ' ' T (^-1) 

W{uo,0, ■ ■ . ,Uo,Ar-l) 

CN,n := [En^I - En){EN-2 - En) . . . {En+1 -En) Tl = 0, . . . , N - 2 

Wn-i{uo.o, . . . ,uo,N-i) VF(wo,o, ■■• ,"0,^-2) „s 

UN N-1 — '■ '■ ^ — — 7 7 ■ l^-^J 

W{uo^o, ■ ■ ■ ,uq^n-i) W(uofi, ■ ■ ■ ,uo^n-i) 
For = 1 the first order transformation operator Li^q is constructed with the 
help of uofi. Therefore applying Q with / = ?Io,o yields 

which obviously agrees with the statement. Applying ^ with / = uo,i we obtain 

W{uofi,uo_i) 

= LifiUo,! = . (A. 3) 

"0,0 

In order to prove the statement for TV = 2 we build the linear SUSY-operator 
-^2,1 from Mi.i and act with it on the unphysical solutions ui,i and ui,o- As first result 
we obtain 

"2,1 = ^2,1"1,1 = = TTT-. r 

"1,1 M/("o,o,"o,i) 
where ([9|) and (|A.3p have been used. For the second function we find 

"2,0 = -^2,l"l,0 = ^2,1 — L'IqUi^i = L+Q_Li.o"o,l 

"0,0 "o,o"i,i ' "o,o"i,i ' 

= {E^-Ea)^^u^,^^{E,-E^)-—^ ^. (A.4) 

"0,0"1,1 VK(uo,o,"o,i) 
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Here the operator £2,1 = ^^^i.iaj/'^i.i+^K has been replaced by Lj^q = —uofix/uo_o — dx 
with the help of (I/uq^o)' = ^^o.o/'"o,o- Moreover, we have used (|A.3|) . the 
factorization rule L^qLi^q — ho — Eq and the Schrodinger equation /io'"o,i ~ EiUq^i. 

Assuming finally that the representations (jA.f p and (jA.2[) are valid for UN,n and 
un,n-i we prove them for ujv+i.„ and un+i,n- 

We start with ujv,jv-i ^Af+LW- To go from hpj to hj^^i only the linear (one- 
step) SUSY transformation Ljv+i.iv is required. Applying ([5]) and combining it with 
the normalization condition W{un^NtUn,n) = 1 and the Crum-Krein formula ((T]) [or 
(|45p] we find the equivalence chain 

_ W^(ujv.JV,WAr^Ar) 1 W {uQfi, . . . , Uq^n -l) 
UN+l.N — ^^N+1,NUN.N — ~ = TT^/ \~ ■ 

un,n unm W (uq^o, . . . ,uo^n) 

Comparison with (|A.2[) shows that the proof is done. 

The proof of the induction u^.n ^ wjv+i.n is less obvious. In section l4?2] the linear 
intertwiners Lm+i,m have been build strictly incrementally from the corresponding 
ground state eigenfunctions Um,m (of the Hamiltonians hm) as transformation 
functions. Here, we need a more general non-incremental construction scheme. In 
order to facilitate it, we first introduce a very detailed notation for general polynomial 
intertwiners L^^q indicating explicitly the energy levels of the transformation functions 
from which they are built. Based on the Crum-Krein formula ([7]) we set 

with Ui e Z+ being any energy level numbers of the Hamiltonian hg. The determinant 
structure of (jA.5[) immediately implies the generalized kernel property 

L'C:,-'"'^uo,n = (A.6) 

and the invariance of the operator ^vj^.^g'"^' "'"''^ with regard to permutations 
(ai, 02, . . . ,ak) ^ cr(ai, 02, . . . , a^) 

^(T(ai,a2,...,afc) _ ^ (ai ,02 , . ...oj.) /. -\ 

Recalling that (|A.5[) can be built from a chain L^kk-i-^k^-i]c-2'''-^'fo' '^^ linear 
intertwiners „ we conclude that the ordering with regard to energy levels is 

inessential in the transformation chain and that we can split it into sub-chains with 
any permuted combination of levels 

a(ai, . . . ,afe) (61, . . . , 6^, ci, . . . , cc) , B + C = k 

r(ai,...,afc) _ J- (ci....,cc) T {hi,.. .,bB) _ r (bi , . . .,6i3 ) t (ci . . ,cc) f \ o\ 

^kfi ^ ^B+C,B ^B,0 ^ ^C+B,C ^C,0 ■ >>^-°7 

Apart from this full commutativity of the transformations we note their associativity 

j-ias) f j(a2) j{ai) \ _ ( j {^1.) t M \ j- (ai) , , 

Commutativity and associativity can be used to re-arrange a sequence of 
transformations in any required order. 

In accordance with the intertwiners L we denote eigenfunctions as 

^(aua,,...,a,) _ ^iai,a,,...,a,) ^^^^ ^ fl, ^ « , * = 1 , . . . , fc . (A.IO) 
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Another ingredient that we need is the representation 

Wjupfi, . . ■ ,Mo,JV-i) _ ( l^N-l-n W{uofi, . . ■ , Up^n-l.UQ.n+l, ■ ■ ■ , Mo,jV--l , UQ.n) 
Wn{uofi, . . . , UQ^N-l) Wn{mfi, . . . , Uq^n ~l) 

( 1 \Ar-l-n 7- (0,...,n-l,n+l,...,Ar-l) 

/ iNW-l-ji (0,...,n-l,Ti+l,...,Af-l) /'A11^ 
= (-1) "7V-1,„ =: VN-l,n (A.ll) 

which immediately foUows from the Crum-Krein formula (IA.5I) and definition (|A.10|) . 
From the physical solution vn-i.u we build the operator^lj 

r(") ._/-(") , 1_ '"jV-l,n,3: ^ 

r(")+ ?-(")+ r 1 Vl\f—i n,x o / A 1 o\ 

and the corresponding Hamiltonian h^^'^^"' with 

r(n)+ 7-(") _ , (0,...,n-l,n+l,...,Ar-l) p / A T?^ 

^N,N--1-'^N,N-1 ~ "-N-l ■ \±^.iO) 

We start the proof of the induction UN,n UN+i,n with the following 
transformations 

_ j-iN) _^ j-{N) Wniuofi, . . . ,Uo^N-l) 

UN+l.n — N ^N.n — N ,n J^N+1 N T777 T 

W [Uo,0, ■ ■ ■ ,Uo,N-l) 

_ t{N) -1 



N,; 



UN,NVN-l,n ^AT-i „ 



= CN,n 'w-l '^N,N ■ (A. 14) 

Un,NVN-1,7i 

In order to obtain the operator product ^v-i^jv jv-i ^^^^ (|45l) . the composition 
rule (jA.Sp and the Crum-Krein formula (jA.Sp 

r(0,l,...,N-l) W{uo.O, ■ ■ ■ ,Uo^n) 

UN.N^-i^MQ "O.JV = 7777 7 

W{uq^q,...,Uo,N-1) 

_ r(") r (0,l....,n~l,ri+l,...,Af-l) 
_ .(n) (O,!,...,"-!,™-!-!,...,^-!) 



where 



(0,l,...,n-l,r! + l,...,Af-l) _ ^ (0,l,...,n-l,n+l,....Af-l) 



(A.16) 



This gives 



/-< 1 r (")+ 

UN.NVN-l.n 



1 



(n)+ T(n) (0,l,...,n-l,n+l,...,A'"-l) 



, „, ^N.N-l^N,N-l^N-l,N 
UN.NVN-l,n 



For C 7^ holds L^^^^_-^ [Cdjv_„] = -i'^^jy— 1 1'"-'^,"] ^° that the sign factor (—1)"'^ " plays no 
role in the operator L,"'' 
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_^ 1 fj^ p ^ (0,l,...,n-l,n+l,-,JV-l) 

^ Wn(Mo,o,...,Mo,Jv) /A 17^ 

= C-Ar+i,„— — (A. 17) 

where the last hne has been obtained by expressing u^^'^'i'j^ i,n+i,...,N i) ^ un^n and 
VN~i,n via (jA.16p . (US]) and (jA.lip in terms of their Wronskian fractions. With (|A.17p 
the proof is complete. 



References 

[1] Groschc C and Steiner F 1998 Handbook of Feynman path integrals (Berlin: Springer) 
[2] Sukumar C V 1985 Supcrsymmctric quantum mechanics of one-dimensional systems J. Phys. 
A: Math. Gen. 18 2917-36 

Cooper F, Khare A A and Sukhatme U 1995 Supersymmetry and quantum mechanics Phys. 
Rep. 251 267-85 

Junker G 1996 Supersymmetry method in quantum and statistical method (Berlin: Springer) 
Bagchi B 2000 Supersymmetry in Quantum and Classical Mechanics (New York: Chapman 
and Hall) 

Mielnik B and Rosas-Ortiz O 2004 Factorization: little or great algorithm? J. Phys. A: Math. 
Gen. 37 10007-35 

[3] Bagrov V G and Samsonov B F 1995 Darboux transformation, factorization and supersymmetry 

in one-dimensional quantum mechanics Theor. Math. Phys. 104 1051—60 

Bagrov V G and Samsonov B F 1997 Darboux transformation of the Schrodinger equation 

Phys. Part. Nucl. 28 374-97 
[4] Samsonov B F and Ovcharov I N 1995 Darboux transformations and non-classical orthogonal 

polynomials Russ. Phys. J. 38/4 58-65 
[5] Aref'eva I, Fernandez D J, Hussin V, Negro J, Nieto L M and Samsonov B F (eds) 2004 Progress 

in Supcrsymmctric Quantum Mechanics special issue of J. Phys. A: Math. Gen. 37 10007—458 
[6] Pershin Yu V and Samsonov B F 2005 Quantum dots created through spherically polarized 

nuclear spins Physica E: Low- Dimensional Systems and N anostructures 28 134-40, 

cond-mat/0401373 

[7] Jauslin H R 1988 Exact propagator and eigenfunctions for multistable models with arbitrarily 

prescribed N lowest eigenvalues, J. Phys. A: Math. Gen. 21, 2337-2350 
[8] Samsonov B F, Sukumar C V and Pupasov A M 2005 SUSY transformation of the Green function 

and a trace formula J. Phys. A: Math. Gen. 38 7557-65, quant-ph/0507160 
[9] Samsonov B F and Pupasov A M 2005 Darboux transformation of the Green's function of a 

regular Sturm-Liouville problem Rus. Phys. J. 48/10 1020—8 
[10] Samsonov B F and Pupasov A M 2006 Exact propagators for complex SUSY partners of real 

potentials Phys. Lett. A 356 210-4, quant-ph/0602218 

Pupasov A M and Samsonov B F 2005 Exact Propagators for Soliton Potentials Symmetry, 

integrability and geometry: methods and applications (SIGMA) 1 20-7, quant-ph/0511238 
[11] Andrianov A A and Cannata F 2004 Nonlinear supersymmetry for spectral design in quantum 

mechanics J. Phys. A: Math. Gen. 37 10297-10321, hep-th/0407077 
[12] Witten E 1981 Dynamical breaking of supersymmetry Nucl. Phys. B 185 513-554 

Witten E 1982 Constraints on supersymmetry breaking Nucl. Phys. B 202 253-316 
[13] Andrianov A A, loffe M V and Spiridonov V P 1993 Higher-Derivative supersymmetry and 

Witten index Phys. Lett A 174 273-279, hep-th/9303005 
[14] Kostyuchenko A G and Sargsyan I S 1979 Distribution of eigenvalues. Selfadjoint ordinary 

differential operators (Moscow: Nauka) 
[15] Levitan B M and Sargsjan I S 1975 Itroduction to spectral theory: selfadjoint ordinary differential 

operators Translations of mathematical monographs; 39 (Rhode Island: AMS) 

Levitan B M 1984 Inverse Sturm-Liouville Problems (Moscow: Nauka) 
[16] Berezin F A and Shubin M A 1991 The Schrodinger equation (Dordrecht: Kluwer) 
[17] Crum M 1955 Associated Sturm-Liouville systems. Quart. J. Math., Ser 2 6 121-127 
[18] Krein M.G. 1957 On a continual analogue of the Christoffel formula from the theory of orthogonal 

polynomials DAN SSSR (Doklady Akademii Nauk SSSR) 113 970-973 
[19] Andrianov A A and Sokolov A V 2006 Factorization of nonlinear supersymmetry in one- 



Exact propagators for SUSY partners 



31 



dimensional Quantum Mechanics I: general classification of reducibility and analysis of tliird 
order algebra Zapiski Nauchnyh Seminarov POMI 335 22-49 
[20] Samsonov B F 1999 New possibilities for supersymmetry breakdown in quantum mechanics 
and second-order irreducible Darboux transformations Phys. Lett. A 263 274-280, 
quant-ph/9904009 

Bagrov V G and Samsonov B F 2002 On irreducible second-order Darboux transformations 

Russ. Phys. J. 45(1) 27-33 
[21] Samsonov B F 2005 SUSY transformations between digonalizable and non-diagonalizable 

Hamiltonians. J. Phys. A: Math. Gen. 38 L397-L403, quant-ph/ 0503075 
[22] Samsonov B F 1996 New features in supersymmetry breakdown in quantum mechanics Mod. 

Phys. Lett A 11 1563-1567, quant-ph/9611012 
[23] Morse P and Feshbach H 1953 Methods of Theoretical Physics, Vol 1 (New York: McGraw-Hill) 
[24] Gelfand I M and Shilov G E 1964 Generalized functions 1-3 (New York: Academic) 

Gelfand I M and Vilenkin N Ya 1964 Generalized functions 4 (New York: Academic) 
[25] Dodonov V V and Man'ko V I 1987 Evolution of multidimensional systems. Magnetic properties 

of ideal gases of charged particles Proceedings of the Physical Lebedev Institute, (Trudy FIAN) 

183 182-286 

[26] Erdolyi A 1953 Higher Trascendental Functions (New York: McGraw-Hill) 

[27] Feynman R P and Hibbs A R 1965 Quantum Mechanics and Path Integrals (New York: McGraw- 
Hill) 

[28] Prudnikov A P, Brychkov Yu A and Marichev O I 1981 Integrals and series. Vol. 1 (Moscow: 
Nauka) 

[29] Abramowitz M and Stegun I A 1964 Handbook of mathematical functions (Washington DC: 

National Bureau of standards) 
[30] Giinther U, Samsonov B F and Stefani F 2007 A globally diagonalizable a'^— dynamo operator, 

SUSY QM and the Dirac equation J. Phys. A: Math. Theor. 40 F169-F176, math-ph/0611036, 



